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Abstract 

In this paper, we try to understand the geometry for a nonlocal nonlinear Schrodinger 
equation (nonlocal NLS) and its discrete version introduced by Ablowitz and Musslimani. We 
show that, under the gauge transformations, the nonlocal focusing NLS and the nonlocal de- 
focusing NLS are, respectively, gauge equivalent to a Heisenberg-like equation and a modified 
Heisenberg-like equation, and their discrete versions are, respectively, gauge equivalent to a 
discrete Heisenberg-like equation and a discrete modified Heisenberg-like equation. From the 
gauge equivalence, although the geometry related to the nonlocal NLS is not very clear, we can 
see that the properties between the nonlocal NLS and its discrete version and NLS and discrete 
NLS have big difference. By constructing the Darboux transformation for discrete nonlocal 
NLS equations including the cases of focusing and defocusing, we derive their discrete soliton 
solutions, which differ from the ones obtained by using the scattering transformation. 
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1 Introduction 


Very recently, Ablowitz and Musslimani investigated a nonlocal NLS equation [T] 

iqtix,t) + qxx{x,t) ± 2q{x,t)q*{-x,t)q{x,t) = 0 , ( 1 ) 

which is derived from a new symmetry reduction of the well-known AKNS system, where q{x,t) is a 
complex valued function of the real variables x and t and * denotes complex conjugation. The nonlocal 
NLS equation ([1]) is a new integrable system possessing the Lax pair, inhnitely many conservation laws 
and it is solvable by using the inverse scattering transformation. Like the classical NLS equation, the 
nonlocal equation m keeps the following parity-time transform invariant 

X—X, t ^—t, q{x,t) ^ q*{x,t). (2) 

Hence, it is PT symmetric [2] and can be regarded as a mathematical model describing wave prop¬ 
agation phenomena in PT symmetric nonlinear media [3 HE]- The nonlocal NLS has attracted 
the attention of researchers since its special properties. For example, by using the inverse scattering 
transformation, Ablowitz and Musslimani obtained its breather solution [I]; In [6] Sarma etal showed 
that the PT-symmetric nonlocal NLS ([T|) can simultaneously support both bright and dark soliton 
solutions; Valchev [7] studied some basic properties of the nonlocal focusing NLS equation including 
its scattering operator, the special solutions by using the dressing method, and the Hamiltonian for¬ 
malism. In addition, dark and antidark soliton interactions in the nonlocal defocusing NLS has been 
discussed [8]. 

On the other hand, Ablowitz and Musslimani also investigated a discrete version of the nonlocal 
NLS O 

W-1" Qn+l + Qn-1 — ‘^Qn QnQ*-TXQn+\ + Qn-l) = 0, (3) 

dT 

which is a discrete PT symmetric model, and it contains a linear Lax pair and an infinite many con¬ 
servation laws. In Ref. 9, a discrete one-soliton solutions with the unique features of power oscillations 
and singularity formation has been derived by using a left-right Riemann-Hilbert formulation. In [6] 
Sarma etal investigated another discrete PT symmetric nonlocal NLS 

i— -h k{an+l + CLn-l) + OnO-n = O' (4) 

dr 

We have known that the focusing NLS and defocusing NLS are, respectively, gauge equivalent to the 
Schrodinger flow of maps from into 5^ in and from R^ into in R^~^^ [in]-[Il]- This gives 
the geometric explanations for the NLS equation. The geometry related to the discrete NLS has also 
been discussed [laiis]. So, what is geometric interpretation for the nonlocal NLS and its discrete 
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version? In this paper, we try to understand the geometry for the nonlocal NLS ([T|) and its discrete 
version ([3|). We will show that, under the gauge transformations, the nonlocal focusing NLS and 
the nonlocal defocusing NLS are, respectively, gauge equivalent to a Heisenberg-like equation and a 
modified Heisenberg-like equation, and their discrete versions are, respectively, gauge equivalent to a 
discrete Heisenberg-like equation and a discrete modified Heisenberg-like equation. From the gauge 
equivalence, although the geometry related to the nonlocal NLS and its discrete version is not very 
clear, we can see that the properties between the nonlocal NLS and its discrete version and NLS and 
discrete NLS have big differences. By constructing the Darboux transformation for discrete nonlocal 
NLS equations including two cases of focusing and defocusing, we also derive their discrete soliton 
solutions, which differ from the ones obtained by using the scattering transformation [9]. 


2 Gauge equivalent structures of the nonlocal focusing NLS and the 
nonlocal defocusing NLS 

In this section, we try to understand the geometry related to the nonlocal focusing NLS and the 
nonlocal defocuing NLS through investigating their gauge equivalent structures. We will show that 
the nonlocal focusing NLS and the nonlocal defocuing NLS are, respectively, gauge equivalent to a 
Heisenberg-like equation and a modified Heisenberg-like equation. From the gauge equivalence, we 
can see that there exist big differences between the nonlocal NLS equation and NLS equation. 

We first recall that the nonlocal focusing NLS equation 

iqt{x,t) -h qxxix,t) 2q{x,t)q*{-x,t)q{x,t) = 0, (5) 

is yielded by the integrability condition of the following linear problem [T]: 

(fx = ipt = Nip (6) 

with 

—i\ q{x,t) ] jY [ —‘^i^‘^ + iq{x,t)q*{—x,t) 2Xq{x,t) + iqx{x,t) 

—q*{—x,t) zA y ’ y —2Xq*{—x,t)+iqx{—x,t) 2iX‘^ — iq{x,t)q*{—x,t) 

Under the following gauge transformation m- 

M = G-^MG-G-^Gx, N = G-^NG-G-^Gt, (7) 

where G is a solution of the system (l6|) for A = 0, i.e., 

Gx = M{0)G, Gt = N(0)G, (8) 
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we can obtain 


M = -iXG-^asG = -iXS, N = -2iX‘^S + XSS^. 


(9) 


The compatibility condition Mt — Nx + [M, A^] = 0 yields a Heisenberg-like equation 

St='-[S,Sxx]. 

From the structure of the matrix M(0) and N{0), we can see that G in ([8|) has the form 


( 10 ) 


G = 


g{x,t) f*{-x,t) 


Hence the structure of the matrix S in equation (|1U|) can be given by 

-1 1 f{x,t)f*{-x,t)+g{x,t)g*{-x,t) 2f*{-x,t)g*{-x,t) 


s = a-‘.,G = ^ 


-‘^f{x,t)g{x,t) 


-f{x,t)f*{-x,t) - g{x,t)g*{-x,t) 


( 11 ) 

where H = f{x,t)f*{—x,t) — g{x,t)g*{—x,t). We should remark here that we say equation (fTO]) is 
a Heisenberg-like equation since it has the same form as the Heisenberg equation, but there exists a 
big difference between (IlOp and the Heisenberg equation. In fact, in the case of focusing NLS, i.e., 
q*{—x,t) = q*{x,t), G in equation ([8|) has the form 


G = 


f{x,t) -g*{x,t) 
g{x,t) f*ix,t) 


and thus S = G ^cjsG has the form 


S = 


1 ( |/|2-|5p -2f*g* 

l/P + bPl -2fg bP-l/P 


Set / = a{x,t) -|- ib{x,t),g = c{x,t) + id{x,t), then the matrix S can be written as 


S = 


Si S2 — ISs 

S2 + iS3 -Si 


where the vector S=(si, S 2 , ss)^ E in R^, and Sj,j = 1, 2,3 is given by 


Si = 


l/P-H' 


2{bd — ac) 

\f? + \ 9 r i/p + isP’ 

So, equation (fTOll reduces to the Heisenberg equation 

Sf = S X Sa;a;. 


S3 = 


—2{bc + ad) 

w+w 


( 12 ) 
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( 13 ) 


For the one-soliton solution of the nonlocal focusing NLS ([5]) 

q{x,t) = csch(8a6t — 2iax)), 


we solve equation (l8|) as follows: 


G = ( Va^+iP 

—b—ia cot]i{Sabt—2iax) 
\/a^+^ 


—6+m coth{Sabt—2iax) 

V 

^^g-26rr+4»(a2-62)tesch(8abt-2ma:) 

VaP+P 


Hence S can be solved as 


•S’—! 


1 + 


2a^ csciP {8abt—2iax) 

aP+P 


2ae^*’“' ^ ^^{—ib+acoth{8abt—2iax))csc]i(8abt—2iax) 

PTiP 


—2ae“^*'^+'^*(“^“*'^)*(ife+acoth(8a6t—2iaa;)) csch(8abt—2iax) 

IPPP 

2a^ csch^ (8abt—2iax) 

aP+iP 


-1 - 


Set / = a{x, t) + ib{x, t),g{x, t) = c(x, t) + id{x, t), then the matrix S in equation (fTTl) can be written 
as 

Si{x,t)+ iS2{x,t) -{S3{-X,t) - iS4{-X,t)) 

y S 3 {x,t) + is 4 {x,t) -{si{x,t) + is 2 ix,t)) 

where 



Si(x,t)si(-a:,t) + S 2 {x,t)s 2 {-x,t) - S 3 {x,t)s 3 {-x,t) - S 4 {x,t)s 4 {-x,t) = 1, 
2si(x,t)s2(x,t) + S 3 {x,t)s 4 {-x,t) - S 4 {x,t)s 3 {-X,t) = 0, (14) 

and Sj{j = 1, 2,3,4) is given by 

si(x,t) = ^, S2(a;,t) = ^, S3{x,t) = ^, S4ix,t) = ^, 

with 

r ={a^{x,t) + b^{x,t)){a^{—x,t) + b^{—x,t)) + {c^{x,t) + (f{x,t)){c^{—x,t) + cf{—x,t)) 

— 2{a{x, t)c{x, t) + b{x, t)d{x, t)){a{—x, t)c{—x, t) + b{—x, t)d{—x, t)) 

— 2{b{x,t)c{x,t) — a{x,t)d{x,t)){b{—x,t)c{—x,t) — a{—x,t)d{—x,t)), 

fill =ia‘^{x,t) + b‘^{x,t)){a^{—x,t) + 6^(—x,t)) — {(?{x,t) + d‘^{x,t)){c^{—x,t) + d^{—x,t)), 
si 2 =2{a{x,t)c{x,t) + b{x,t)d{x,t)){b{—x,t)c{—x,t) — a{—x,t)d{—x,t)) 

— 2{b{x, t)c{x, t) — a{x, t)d{x, t)){b{—x, t)d{—x, t) + a(—x, t)c(—x, t)), 

5 21 =2(a^(x, t) + b^(x, t))(d(x, t)b(—x, t) — c(x, t)a(—x, t)) 

+ 2(c^(x, t) + d^{x, t)){a{x, t)c{—x, t) — b{x, t)d{—x, t)), 

5 22 = — 2(a^(x, t) + b^{x, t)){d{x, t)a{—x, t) + c(x, t)b{—x, t)) 

+ 2{c^{x, t) + d'^{x, t)){b{x, t)c{—x, t) + a{x, t)d{—x, t)). 
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We thus see that equation (fTOt) can be transformed into 


— 2 ('®4'S3xx( 3;) S3( x'jS^xx 'S3xx^’4( 'S3'®4xx( 3^)) , 

—2 (®3®3a:a;( x') S3( X^S^xx “b S4®4xa;( x') S4( x)s4xx) 1 

dss 

j, — S 2 S 3 XX S 3 S 2 XX “b 'Sl'S4xx ^’45111) 

at 

dS4 

], —S3'Slxx ®lS3a:a; “b S 2 S 4 XX S4S2xx- 

at 


( 15 ) 


Next we will discuss the gauge equivalence for the nonlocal defocusing NLS equation 

iqt{x,t) + qxx{x,t) - 2q{x,t)q*{-x,t)q{x,t) = 0 , (16) 

which has the Lax pair [1] 


(fix = M^p, cpt = Nip 


(17) 


where 

M = 


A q*{-x,t) 
q{x,t) -A 


N = i 


—2A^ + q{x, t)q*{—x, t) —2Xq*{—x, t) + q%{—x^ t) 
—2Xq{x, t) + qxix, t) 2A^ — q{x, t)q*{—x, t) 


Under gauge transformation: 

M = G-^MG-G-^Gx, N = G-^NG-G-^Gt, 

where G satisfies 


Gx = M(0)G, Gt = NiO)G, 


we can obtain 


M = AG-V 3 G = iXS, N = 2X^S + iXSSx 


(18) 

(19) 

( 20 ) 


with S being defined by 5 = —iG ^cr^G. The compatibility condition Mt — Nx + [M, iV] = 0 yields a 
modified Heisenberg-like equation 


St = ^[S,Sxx]. 

From the structure of the matrix M(0) and A^(0), we can see that G in (|19p has the form 


( 21 ) 


G = 


f{x,t) -g*{-x,t) 
9{x,t) f*{-x,t) 
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Hence the matrix S is given by 

S = -( Six,t)g*i-x,t) - f{x,t)f*{-x,t) 2f*{-x,t)g*(-x,t) 

2 /(x,t) 5 f(x,t) f{x,t)f*{-x,t)-g{x,t)g*{-x,t) 

where H = f{x,t)f*{—x,t) + g{x,t)g*{—x,t). We remark here that we say equation ([ 2 T]) is a modified 
Heisenberg-like equation since it has the same form as the modified Heisenberg equation, but it has 
big difference with the modified Heisenberg equation. In fact, we see that in the case of defocusing 
NLS, the matrix G in equation (|19p has the form 

fix,t) gix,t) \ 

V -g*{x,t) J ’ 

and thus S = —iG~^a 3 G has the form 

i [ -ilmjfg*) -\g\‘^ \ 

Re{fg*) y -|/|2 ilm{fg*) ) ' 

Set / = a{x,t) -b ib{x,t),g = c{x,t) + id{x,t), then the matrix S can be rewritten as 



S = 


Si i{s3 - S 2 ) 

*(S2 + S3) -Si 


where the vector S=(si, S 2 , ss)^ G in i.e., s\ + s\ — s^ = —1, and Sj,j = 1, 2, 3 is given by 


be — ad 


Si = 


S2 = 


|9P-|/1^ 


S 3 = 


-l/P-H' 


ac + bd^ ^ 2{ac + bd)' 2{ac + bd) 

Thus, the modified Heisenberg-like equation (1211) leads to the modified Heisenberg equation 


St = 


SxS, 


(23) 


where x denotes the pseudo cross product in defined by axb = (0263 — 0362,0361 — 0163 , —(0162 — 

0261 )). 

For a solution q{x, t) of nonlocal defocusing NLS given by 

q{x, t) = sech( 8 a 6 t + 2i6x)), (24) 


the solution to equation (1191) is 


G 


—a—ibtanh{8abt+2ibx) 4,(0^ sech(8abt+2ifex) 

—sech(8abt+2ifea;) —a+ibtanh{8abt+2ibx) 

\/a'^+b‘^ 
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and 


1 2b^ sech^ (Sabt-^2ibx) 

5=1 22 

2be^<°- sech{8abt+2ibx){a+ibta.nh{8abt+2ibx)) 

aF+IP 


2j,e-4i(a^-6^)t+2ar^ sech{8abt+2ibx){a-ib tanh{8abt+2ibx)) 
a^+fe^ 


-1 + 


2b^ sech^ (8abt+2ibx) 

aF+iP 


Set / = a{x,t) + ib{x,t), g{x,t) = c{x,t) +id{x,t), then the matrix 5 in equation ([22]) can be written 


as 


where 


I Si{x,t) + iS 2 {x,t) -{S 3 {-X,t) - iS 4 {-X,t)) 
y S 3 {x,t) + iS 4 {x,t) -{si{x,t) + iS 2 ix,t)) 


Si{x,t)si{-x,t) + S2{x,t)s2{-x,t) + S3{x,t)s3{-X,t) + Si{x,t)s4{—X,t) = 1, 
2 si(x,t)s 2 (x,t) + S3{x,t)s4{-x,t) - S4{x,t)s3{-X,t) = 0 , ( 25 ) 

and Sj{j = 1, 2,3,4) is given by 

si(x,t) = ^, S2ix,t) = ^, S3{x,t) = ^, Si{x,t) = ^, 

with 


r ={a?{x, t) + t))(a^(—X, t) + 6 ^(—x, t)) + (c^(x, t) + d?{x, t)){c^{—x, t) + cf{—x, t)) 

+ 2 (a(x, t)c{x, t) + 6 (x, t)d{x, t)){a{—x, t)c{—x, t) + b{—x, t)d{—x, t)) 

+ 2 ( 6 (x, t)c{x, t) — a(x, t)d{x, t)){b{—x, t)c{—x, t) — a{—x, t)d{—x, t)), 
sii = 2 (o(x, t)c{x, t) + b{x, t)d{x, t))(a(—x, t)d{—x, t) — b{—x, t)c{—x, t)) 

+ 2{b{x,t)c{x,t) — a{x,t)d{x,t)){b{—x,t)d{—x,t) + a{—x,t)c{—x,t)), 
si 2 = - (a^(x,t) + b‘^{x,t)){o?{-x,t) + 6 ^(-x,t)) + {c^{x,t) + d‘^{x,t)){c^{-x,t) + d‘^{-x,t)), 

5 21 = — 2 (a^(x, t) + 6 ^(x, t)){d{x, t)a{—x, t) + c(x, t)b{—x, t)) 

— 2(c^(x, t) + d^{x, t)){b{x, t)c{—x, t) + a(x, t)d{—x, t)), 

5 22 = 2 (a^(x, t) + 6 ^(x, t))(c(x, t)a{—x, t) — d{x, t)b{—x, t)) 

+ 2 (c^(x, t) + d^{x, t)){a{x, t)c{—x, t) — b{x, t)d{—x, t)). 


Thus, equation ^T\i can be transformed into 


—2 ('53'S3a;a;( x') S3( x')S3xx T S4S4.xx(^ S4( x)s4xx) j 

(^4S3xx( x) S 3 ( x)s4xx T 'S3a:3;'S4( x) 'S3S4a;3;( x)) , 

dS3 

1 , — ^l^Sxx “T S 2 S 4 XX ^4^2xx'f 

at 

dS4 

—7— =S4Sixx — •SlS4a;a; + ■S 3 S 2 XX ~ S2S3xx- 

dt 


( 26 ) 






In summary, although the geometry related to the nonlocal focusing NLS and the nonlocal defocusing 
NLS is not very clear, we can see, from their gauge equivalence, that the properties between the 
nonlocal NLS equation and NLS equation have big differences. 


3 The soliton of discrete nonlocal NLS and gauge equivalence 


In this section, we will seek the soliton solution of the discrete nonlocal focusing NLS and the discrete 
nonlocal defocusing NLS through constructing their Darboux transformations. We will show that 
there is no singular point in the discrete one-soliton solution, which is distinguished from that given 
in [9]. We will also show that the nonlocal discrete focuing NLS and the nonlocal discrete defocuing 
NLS are, respectively, gauge equivalent to a discrete Heisenberg-like equation and a discrete modified 
Heisenberg-like equation. From the gauge equivalence, we can see that there exist big differences 
between the nonlocal discrete NLS equation and discrete NLS equation. 

3.1 The soliton of discrete nonlocal focusing NLS and gauge equivalence 
The discrete nonlocal focusing NLS is as follows: 

i-^ - 1 “ Qn+l + Qn-l — ‘^Qn + QnQ*-n{Qn+l + Qn-l) = 0, (27) 

cLt 

which has the discrete Lax pair 

Tn,r ( 2 ^) 

with 



2 ; QLnZ ^ \ 

-QnZ Z-^ j ’ 

1 - -h z - Z~^ - Qf^Qn-1 
Qn T Qn—lZ 


-Qtn + Q*-n-iZ-^ 

— 1 + 2 + Z — Z QnQ-n-1 


We remark here that the Lax pair (|28p is different from the one given in [9]. Equation (|27l) is PT 
symmetric similar to the classical integrable discrete NLS. Introduce (pn^ = where the matrix 


with constraint condition: 


r„ = 


2 ; -I- anZ ^ bnZ ^ 
CnZ dnZ -I- Z~^ 


— C 


* 

— 725 




d 


* 

— 72 * 


(29) 


(30) 


Suppose Qn is a solution of (ITH) and tpn = (v^i.n, V^ 2 ,n)^ is an eigenfunction of linear problem (ESP with 
z = zi- Then one can check that V’n = {—T 2 -mTi -nY' eigenfunction when z = (z]")”^. 
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Assume that detrn,( 2 ;i) = 0, then the two column vectors in Tn{zi){ipn,'tpn) are linear dependent. 
Thus we get 


Cn — 


zf - ^)Tn 

1 + TnT* 


dn, — 


'Tn'T—n 
1 + Tnr*„ 


where = ^P 2 ,nl^i,n- we have proved that the new linear problem 

where 

MW = Tn+lMnT-\ iVW = ^Tn,r + T„iV„)r-l, 


(31) 


(32) 

has the same form as the linear eigenfunction equation (f28]l except that the Qn,Q*-n ia Mn,Nn are 
replaced by Q^n , Q*-} in . 

The relation between old potential Qn and new potential Q^n is 

QW = Qndn+l - Cn+l- (33) 

For the seed Qn = 0 and zi = a + i/3, the eigenfunctions are (fi^n = Zie^'^,ip 2 ,n = with 

^ = i(l - ^^1 + - z^^),r] = i(-l + z^'^ + zi - z^^). 

So the discrete one-soliton solution is 

^_2(n+l)(^_2 _ 


qw = 


Its norm is 


_(a2 + / 32 )-n|(g ^2 ^ ^2)2 _ i|g2«/3((a^+/3^)-^-l)r_ 

((a^ + /3^)^ + e4Q:/3((a2+/32)-2_i).,- sin^ u 

where v = 2(n — l)Arg(a — i/3) — 2(n + l)Arg(Q; + i/3). 

If the zi is the form of polar coordinates, i.e., zi = ae*^ and a / l,/3 E (—7r,7r], then we get 




[ 1 ] 


(34) 


(35) 




a — lie 


- a —T ]T 


”2—^ '-sin(2^) 


_4( 0,2-4 


-A'It sin(2/3) , —2(q 2—h')Tsin(2/3) . 

+2Q;'^e V cos(4/3n) 


(36) 


a“ + e 

Note that there is no singular point in discrete one-soliton (I34p . which is distinguished from that given 
in [^. Fig. 1 gives the shape of the discrete one-soliton solution with a = y/5j2,l3 = arctanl/2. 

On the other hand, if we take the seed solution Qn = p,(f) ^ R, then solving linear 

isospectral equation (|28]) yields the following eigenfunctions: 

Pin = ip{z-^ - l)C”e^i" + (A2 - q)D^e^^\ 

(37) 

(P2„ = (c-e^i"(Ai -p) + ip{z^ - l)74”e^2") , 
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where 


Ai — 


A 2 — 


C = 


i{l -z‘^ + 2z^ - 2p^z^ - 2z) + (z^ - 1)\/A 
2 i 2 ’ 

i(l -z^ + 2z^ - 2p^z^ - 2z) - (z^ - 1)\/A 
2 i 2 ’ 

z(z^ - 1 - 2 p 2 - iVA) 


z^ + 1 + i\/A 


2z 


D = 


— 1 — i\/A 


A = Ap^z^ - (z^ - 1)^, 


p = i(l-z^ + z-z ^ - p^), q = i(z ^-h z - z ^-1-p^). 
So, new soliton solution is 


( 38 ) 


n[l] = (^1 ^ zf )Tn+l _ 

1 + r„+ir*„+i 




1 + Tn+lT'^n+l 


(39) 


where 


with 


^n +1 — 


-pZ^e'^^P + Tn 

Z^ + pe~'^^P^'^~'^^Tn 


_ ^2ip^r+iP - 1) + \/A + 2z/3Z^g"eA2 ^ 1 )^ 
2ip + fi(z2 - 1) + \/a) 6'’^eA2-^i)^ ’ 


Tn = -e 


.V-2p^-l-iVS) ^ - A. = - 1)VA, 

C z‘^{l + 2p‘^) — 1 — iy/A 


A2 - g = -(Ai -p) = 


(1 - z 2 ) (^i (^2 _ 1 ) + ^ 


2^2 


It is interesting to note that by choosing the proper parameters, q'^u can be rewritten as 


(40) 


(41) 


QW=pe2*^'"+*^/(n + 7r,n-/3T) 


(42) 


This means that the solution \Q^n\ is a traveling soliton solution. In the specihc case, we give out the 
shapes of |(5n^ |. Set 2 : = a + i6, + 6^ / 0, / 1, then we get A = — 1 — — ^2(2 + 52 + 4p2) -\-a^[2 + 

66 ^ + 4p2) + 4a6i(l — + 6 ^ + 2 p 2 ). In the case of ImA = 0, i.e, a = 0or6 = 0orl — + 6 ^ + 2 p 2 = 0 , 

we discuss the properties of Qn^. 

Case 1: 6 = 0 

• When ReA < 0, i.e., 2 |p 2 ;| < \z‘^ — 1|, we have 



_g*(2p^T+(/>) 


,2(z- 


e2{z 2 


+ e(^ 2 
i)V^iTM4^ + e(^“^ 


i)^^"M 2 n + M 3 
+ Me ’ 


(43) 
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where 


Min =2e^^+^p^z‘^{l + z^){l -z^ + \/^), 

M2n = - pM^n = 2e^pz\l - z^) - 0^ + (1 + 0^)./^ + (1 + 2p‘^){z‘^0^ - 1)) , 

Ms =20(1 + (1 - + 3(1 + p‘^)z‘^ - (3 + P^)z^ + [p^z^ - (z^ - 1)^) V-A , 

M 4 n =9 ^^Mq, Mg = 2p9z‘^{l + z^){z^ - 1 + \/^). 

Note that Mjn,j = 1, 2,..., 6 is independent of r, so |(5n^('r)| is a discrete time-period soliton solution 
with T = (^-2 • The discrete soliton solution \Qin{T)\ is shown in fig. 2 where the parameters 

z = 2, p = 1/2, which is a nonzero background solution with period T = To describe the back¬ 
ground clearly, we give the graphs at three different space lattice points n = —5, n = 0 and n = 5, 
respectively. 


When A = Re/S. > 0, we can get a breather-like solution on space lattice n 


QW(r) = _eh2pV+0) 


e2(^-"-i)^^M4n + + MQn ’ 


where 


Min =Z^M^n = ‘2pz^{l + z‘^)0^ (^2(1 + p2);j2 _ ^ ^ 

M2 = - pMs = 2p'^z\l - z^) ((1 + ^2 - iVA)0* + (1 + ^2 + iVA)0) , 

M4n =Z^M^n = V^^(l + Z^)0^. 

Fig.3 shows the breather-like solution where z = 2, p = 1. 

Case 2: For 1 — -|- 6^ -|- 2p^ = 0 

In this case, we have Re/S > 0. When 6 / 0, we find that |(5n^(r)| is a two-solitons without singular 
and it has local maximum value (see hg.4, where a = 2,h = p = 1). 

Case 3: For a = 0, he., z = ib 

In this case, one can check Re A < 0. From the formula of we see that the IQn^l is also a discrete 
time-period soliton. Fig. 5 gives the discrete period soliton with T = where b = —1/2,/j = 1. 
Next we will discuss the gauge equivalence for the discrete nonlocal focusing NFS equation. Let 
Sn — Gn^a^Gn, where Gn satisfies the linear problem 


Gn+l = Mn{l)Gn, Gn,r = A,(1)G, 


with the form of 


Gn = 


fn -g*-r 

9n fSn 
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Under discrete gauge transformation 


Mn = G;Xl^nGn, Nn = G-^N^Gn - G^^Gn^r, 


(45) 


we obtain 


Mn =G-^M-\l)MnGn = 




Nn=G-\Nn-Nn{l))Gr^ 


I ( I-z‘^ + z-z^ (z2 - 
" I {z^-l)Qn-l -1 + Z-2 + Z-Z-1 


G. 


=i{z — z ^)I + i 


-iw , .(z-^ + z-^ 


2 ^ 


I ^ Q-l I -Q-n-1 


2 

„-i 


=i{z — z ^)/ + U 1 — 


Here we have used the identities 


Qn—l 1 
+ 2;- 


Qn—1 1 
Gn 


Gr, 


Sn + Sn-1 .Z^ - Z ^ I + Sn-lSn 

— l- 


1 + itr(5n5'„_i) 2 1 + itr(5n5'„_i) 


l+2^r(5"+i'^")“ 1 + Q„gi„’ 


G 


-1 


1 Q-n-l 


Qn— 


^—n—l I /~i _ /^—1 ^ 

I ~ 1 1 1 

1 i / i + n 


+ Sn-lSn 


tr(SnSn-l) 


G 


-1 


1 -Q-n-l 

Qn—l 1 


Gn = 


Sn + 'S'n-1 




Then by using the discrete zero curvature equation = ^n+iMn — MnNn and comparing the 

power of z, we get a discrete Heisenberg-like model 

dSn 


Jn . [*5,1+1,-S', 

= I 


n\ [<5,1, * 5 , 1 - 1 ] 


dr 1 + itr(5n+i5n) 1 + ^tr(5„S'n-i) 
where the matrix Sn is given by 

1 


Sn — Gn Cr^Gn — 


fnf-n + 9n9-n y 
Set fn = cin + ibn, gn = Cn + idn, then Sn has the form 


fnfln - 9ng-n -V-n9-n 

-2/n5n 9ng-n “ /n/-. 


(46) 


( 47 ) 


* 5 „, = 


_ / Sin + iS2n 'S3(-n) “ ^•S4(-„) 


S3n~\~iS/ln ~ {sin iS2n) 
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where 


SinSl(_n) + S2nS2{-n) + S3n'S3(_n) + 'S4nS4(-n) — 1) 
2sin'S2n “1“ 'S3(— n)®4n — 0, 


( 48 ) 


and SjnU = 1, 2, 3,4) is given by 


ai 


-.12 


c21 


C.22 




^2n — 


^3n — 


— 


with 


Tn —{(^n + + ^-n) + i^n + + ^-n) + 2(art?^-n — a-nbn){Cnd-n ~ C-ndn) 

+ 2((Z72(Z_jt, -|- 67 ^, 6 —^^)(c72C_-n, H“ n)? 

={o-n + ^n)(®-n + ^-n) ~ (^n + dn){C-n + 

Sn — ‘^{finO'—n 4“ bnb—n}{Cnd—n C—ndn) 4“ 2(C}^C_}^ 4“ dnd—ji){ci—nbn ®n^—n); 

Sjj — ‘2{p,n 4“ n C-nO^—n) 4“ 2(Cjj 4“ dj^)(&nrf—n ^nC—n) 

Sjj = 2(cijj 4“ b^iOi—ndn 4“ b—nCn) 2(c^ 4“ d-i^iOind—n 4“ bnC—n). 


The matrix equation (j46l) can be rewritten as 



■^ndd'n 4“ B-i-i^F-n 

dnd-n DnJn 

dS2n ^nd^n 

CnJyi dDrilyi 

dr 

Ain 

A2n 

dr Ain 

A2n 

ds^n 

■^nGn 4“ 

Cfid^n “1“ dDfiLyi 

-^ndin 

C'nd'n dDnddn 

dr 

Ain 

A2n 

dr Ain 

A2n 


(49) 
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where 


Ain — ^n + -®n) ^2n — C'n + ^n'l 

= 1 + SlnSl(n+l) ~ 'S2nS2(n+l) + 2 ('®3n'S3(_n+l) + 'S3(-n)S3(n+1) + 'S4n'S4(_n+l) + 'S4(-n)S4(n+1))5 
Bn = S2n'Sl(n+l) + SlnS2(n+l) “ 2 ('®3»^'®4(-n+l) + 'S4(-n)S3(n+1) “ 'S4nS3(_n+l) ” S3(_„)S4(„+i)); 

Cn = 1 + Sln'Sl(n-l) “ S2nS2(n-l) + 2 ('®3(-n)'S3(n-1) + ^SnSs^-n-l) + 'S4(n-1)'S4(-n) + 'S4n'S4(-n-l))5 

1 

Bn 'Sln'S2(n—1) “1“ ®2n'Sl(n—1) 2 (^3n'S4(—n—1) '53(—n)S4(n—1) “1“ '53(n—1)'S4(—n) '54nS3(—n—1))i 

En = S3n'S4(-n+l) “ 'S3(n+1)'S4(-n) “ S3(-n+l)'S4n + 'S 3 (_„) S 4 („+i); 

En = S3„S3(_„+i) — S3(n+l)S3(_n) + S4(-n+l)'S4n “ 'S4(_„)S4(„+i); 

Gn = 2(s2(n+l)'S3n “ 'S2nS3(n+l) + 'Sl(n+l)'S4n “ 'SlnS4(n+l))! 

Hn = 2(sinS3(n,+l) “ 'Sl(n+l)'S3n + 'S2(n+l)'S4n ~ 'S2nS4(n+l))i 
In = —S3nS4(-n-l) “ S3(-_„)S4(„_i) + S3(n-1)S4(-n) + 'S3(_n-1)S4n! 

In = 'S3(-n)'S3(n-l) “ 'S3(n)'S3(-n-l) + 'S4(n-1)'S4(-n) “ 'S4(-n-l)'S4n; 

Kn = 2(s2nS3(n-l) “ 'S2(n-l)'S3n + 'SlnS4(n-l) “ 'Sl(n-l)■S4n); 

En = 2(si(„_i)S 3„, — 'SlnS3(n-l) + 'S2nS4(n-l) ” S2(n-1)'S4n)• 

(50) 

We remark here that if set /*„ = fn,g-n — 9n fn = o,nE ibn^dn = Cn + idn, then we can show 
that the matrix Sn possesses the form 


Sn 


Sin S2n '^S^n 

S2n “1“ is^n 


(51) 


where the vector Sn={sin, S 2 n, S 3 n)'^ £ <5^ in R^, and Sjn,j = 1,2,3 is given by 


Sin — 


l/nP - bn 


l/nP + bnP’ 

and equation (H6l) reduces to 


S2n = 


2(^bndn OinC-n) 

l/nP + bnP 


S2tn — 


2(bnCn “h dndn^ 

l/nP + bnP 


s„, = 2 


5n+l 


X Sr 


1 + S„4-1 • S 


- 2 


Sfi X Sfi_i 

1 + Sn • Sn-1 ' 


(52) 


This means that the equation (|46n is a discrete Heisenberg equation under the condition = 
fni9-n — 9n which agrees with a well-known fact that the discrete focusing NTS is gauge equivalent 
to the discrete Heisenberg equation m- 
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3.2 The soliton of discrete nonlocal defocusing NLS and gauge equivalence 
For the discrete nonlocal defocusing NLS 


its Lax pair is: 


i-^ - 1 “ Qn+l + Qn-l — 2,Qn — QnQ*-n{Qn+l + Qn-l) — 0 , 

cLt 


Tn,T — 


(53) 

(54) 


with 


* ^-1 


Z Q-n^ 
\ QnZ Z-^ 


Nn =i 


1 — Z +Z — Z + Qf^Qn-1 —Q-nTQ*-n-l^ 


-2 


Qn Qn—lZ 


— l-\-Z +Z — Z — QnQ-n-1 


We also remark here that the Lax pair (1541) is different from the one given in [^. We have obtained 
the following Darboux transformation for discrete nonlocal defocusing NLS equation (|53l) : 




with constraint condition: 


Further, we have 


Z + UnZ ^ bnZ ^ 
CnZ dnZ + Z 


bn — — ^—r 


-1 


Tni 


(55) 


(56) 


Cr7, — 


zf - Z;^ ^)Tn 
1 - TnTl^ 


dn — 


'’'nT_n 

1 - TnT*_^ 


and the relation between old potential Qn and new potential Qn^ can be written as 

Qn^ = W+l + Qndn+1- (57) 

For seed solution Qn = 0, one can get discrete 1-soliton solution for discrete nonlocal defocusing NLS 

Z-2{n+l)^zf - 


qw = 


(58) 


1 _ 

where ^ = i{l — zf + zi — zf^),ri = z(—1 + zf'^ + zi — zf^). Set zi = and a 7 ^ 1, /3 E (—tt, tt], then 
we get 

«-2n|^4_i|g-(«^-^)-in(2/3) 






( 59 ) 
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Its shape is given in Fig.6 with a = y/5j2,l3 = arctanl/2. The solution (|58l) differs from one of the 
classical discrete defocusing NTS. Taking the seed solution 


and solving equation 


Qn = pe p,(l)eR 


give the following eigenfunctions: 

= ipiz-^ - + (A2 - 

<^2n = -p) + ipil - ^ 


where 


Ai — 


A2 — 


C = 


i{l -z^ + 2z^ - 2p^z^ - 2z) + (z^ - 1) VA 

2i2 ’ 

i(l — z^ + 2z^ — 2p^z^ — 2z) — {z^ — 1)-\/A 
2i2 ’ 

z{z^ - 1 + 2p2 - iyfK) 


1 + iVA 


2z 


D = 


Z‘ 


— 1 — i\/A 


A = - Vz" - {z^ - 1)^ 


p = i{l-z^ + z-z ^+p^), q = i{z “^ + z - z ^-l + p^) 
We thus obtain new soliton solution 

,*2 .^-2^_ . _^*2 X ^-2 


[1] _ {zf - z^ ^)Tn+i __ 2 ip 2 ^+i 4 > -zl + -r^+ir*n+i 
^ ^ 'T pc -I ^ ’ 

1 7n+ir_„_l_]^ 1 Tn+lT_^_^_i 


with 


"^71+1 — ^ 


-2ipV+i<A - 1) + VA - 2ipz2pri+lg(A2 Ai)7 


2ip + (i{z^ - 1) + VA ) d^+^e(^2-M)r 




D z'^{z‘^+ 2p‘^ - 1 - i'/A) 2 -,\ fX 

77 = -77 WT ^ \2-Xi = {z ^ - 1)VA, 


’(1 - 2p2) _ 1 _ 


h-q = -(Ai -p) = 
)[!] 


(1 — (i{z^ — 1) + a/A 

2^2 


(60) 


(61) 


(62) 


(63) 


Let us discuss the properties of Qh . Set z = o + i6, a2 + ^2 ^ 0, / 1, then we have A = —1 — m — 
6^(2 + 52 — 4p2) + a‘^(^2 + 6^2 — 4p2) _|_ 4abi{l — + b'^ — 2p^). Suppose ImA = 0, i.e, a = 0 or 6 = 0 

or 1 — o2 + ^2 _ 2p2 = 0. 

Case 1: For 6 = 0 

In this case, due to Re A < 0, we find that |Qn^(r)| is a even discrete soliton on time r when 
|(5n^(r)| does not exist singularity for some parameters. Fig.7 gives the discrete soliton |(5n^(r)| where 
b = 2,p = —0.5,4> = 0. Note that IQn^T")! is related to </>. When </> / 0, |Qn^(r)| has a shift along r. 


_ „.-.2pV+i^ e2(^-^-i)v^--Mi^(r) + e(-~^-i)v^--M2.(r) + Ms^^jr) 


QIiKt) = pe 




(64) 
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where 


Min{T) =2/5z202n+l ^ _^2)j _ 

M 2 „(r) =r [(2^(1 - e‘^z^) + 2l2(l - e'^z^) + - l)) /^ + + V) 

M 3 „(r) =y40(A + - 2z‘^p^ - z^p^), 

Min =e^'^Me = 2Ae^+‘^^z‘^{l + z‘^)p‘^, 

M^n =pz^O^ {A\9^ - 1) - V(/02 _ ^ 


with A = z'^ — 1 + ^/—A. 

Case 2: For a = 0 

When ReA > 0, by choosing proper parameters, we get a discrete breather-like soliton (see fig.8) 



= g-h2p^T+(?i) 


g2{z-^ 


i)VAtJ r e(^ 2 

_g g(^-2 


+ Man 


(65) 


where 

Min =pz‘^9^ [‘iz'^p'^ — {Az‘^ — 2i{z‘^ — 1)) A*] , 

Man = {2iAz‘^p^ + {^- - 2iz*p‘^)A*) 9* + 2z‘^p^9{-2 + iAz^ + 2z^ - iA*), 

Msn = 90 *"" {2iA{z^ - 1) - 4zV^ + AA*) ; 

Min =9^z^{AA* - 4zV^), Men = z‘^9*{Ap^ - AA*), 

M^n =2iz^p ((al* - Az^)9* + 9{z'^9* - A)) , 

with A = i{z‘^ — 1) -|- \/A. 

Case 3: A 7^ 0. In this case, we can get a discrete 2-soliton solution with nonlocal maximum value 
when a = 0.5, 6 = 1, and p = 2 (see fig. 9). 

Next, we will discuss gauge equivalence of discrete nonlocal defocusing NLS. Set Sn — —iGn^a^Gn, 
where G„ satisfies the linear problem 


Gn+l = Mn(l)Gn, Gn,r = Nn{l)Gn 


with the form of 



Under discrete gauge transformation 


Mn = G-l^MnG 


Nn = G-^NnGn - 


( 66 ) 
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we obtain 


z + z .z — z ^ 


Mn = - iG-^M-^il)MnGn = - 1 + i 

Nn =G-HNn - Nn{l))Gn 




^ 1 / l-Z^ + Z-Z ^ {Z ^-l)Q*_n-l 

' " I {l-Z^)Qn-l -1+z-^ + Z-Z-^ 


Gr. 


=i{z — z ^)I + i 


-iw . .(z^ + z-^ 


2 ^ 


.2 _ .-2 
+ i—^G-^ 


2 

-l^ 


=i{z-z-^)I + 


-1 -Q-n-l 

Qn—\ 1 

Z ^ — 2 Sn + Sn— 


Qn—\ 1 

\ G„. 


Gr, 


n “1“ *^n—1 .2^ Z I Sn—lSn 

— l- 


1 - itr(5„5n-i) 


2 1 - ltr{SnSn-l) 


Here we have used the identities 


1 


G, 


G. 


-1 I 1 Q-n-l \ 
^ Qn-1 1 j 

-1 / 1 -Q-n-l 

\ Qn-1 —1 


n — ^n-l^n — “ f 


/ - Sn-lSn 


Gr). — 


1 - 5tr(S'„5n-i) 
Sn + Sn-1 


1 - itr(5n5n-l) 


Then by using the discrete zero curvature equation Mn^r = -^n+i^fn — M^Nn and comparing the 
power of z, we get a discrete modified Heisenberg-like model 


dt 




1 - kr(5n5n-l) 


[5’n+l,5n] 

1 - itr(5n+i5n)’ 


where the matrix Sn is given by 


Sn — —iGn^O'sGn — 


fnQ-n + 9nfZ 


9nfln - fn9-r 

-2fnr-n 


-25n5-n 
fn9-n - 9nf-n 


(67) 


( 68 ) 


Let us further investigate the structure of the matrix Sn- Set fn = an + ibn,9n = Cn + idn, then Sn 
has the form 


Sn 


■Sin + iS2n S^n + is^n 
S3n + iS4n —{sin + iS2n) 


19 



where 


'Sln'Sl(—n) “1“ ^2nS2(—n) 'S4nS6n — Ij 

2sin'S2n “1“ 'SSn^Gn “1“ — 0; 


and SjniJ = 1,2, 3,4, 5,6) is given by 


Sin — 


Sin — 


ai 


-.22 


S2n = 


SSn — 


,12 


„31 


S3n = 


Sen — 


c21 


r ’ 

^32 


with 


Tn —(fln + bn){c^n + '^-n) + (®-n + ^-n)(Cn + '^n) + 2(a„6_„ + a_„6„) + C^ndn) 

Sfi — “^{bni—n O^nOi—n^ic-nd—n 4“ C.—ndn) “^{CnC—n dnd—n){Qi—nbn 4“ 0,nb—n\ 

s}? =(a-n 4- b^-n){Cn + dl) - + d^_n), 

Sfi —^(®n 4“ b^) (b—nC—n (^—nd—n) 4“ 2(q,_jj + b_^'){Ojndn bnCn)i 

Syi ~2(Ojj 4“ by^{jl—nC—n 4“ b—nd—n') 4“ 2(fl_jj + b_y^{jlnCn ~\~ bndn), 

Syi ~‘^i.('n 4“ d^{jl—nd—n b—nC—n) 4“ 2(c_jj 4“ d_^(bnCn 0,ndn\ 

Sn^ + d^){a-nC-n 4" b-nd-n) 4" 2(c?_„ 4" d‘^^){anCn 4" bndn)- 


ds\n 

Cfiln “1“ dJfiJ'fi 

^nEn 4- BnFn 

dS2n _ 

J D T 

^n-^-n 

AnFn — BnFn 

dr 

Ain 

A2n 

dr 

Ain 

l> 

to 

ds^n 

Cnd^n 4“ Dndin 

^nPn 4- BnHn 

dS^^yi 

Cnd'n 

-4ni^n - -Bni^n 

dr 

Ain 

1> 

to 

dr 

Ain 

l> 

to 

ds^n 

_ CnXn 4- DnYn 

AnRn 4- BnWn 

dSQn 

_ CnYn — DnXn 

^nld^n BnRn 

dr 

Ain 

A2n 

dr 

Ain 

1> 

to 


( 69 ) 
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with 


^In = C'n + -^n) ^2n = ^ 

=2 (2 “ ■Sln('Sl(n+l) “ Sl(n-l)) + 'S2n('S2(n+1) “ S2(n-1)) ~ S3(n-l)S5n 
—'S3„S5(„+i) + S4(n-l)S6n + S4n'S6(n+l)) ! 

Bn =-^('Sln'S2(n-l) + S2n'Sl(n-l) + S2nSl(n+l) + SlnS2(n+l) + 'S5nS4(n-l) 

+ S4nS5(n+l) + -SenSsin-l) + 'S3nS6(n+l))! 

Cn ~ + 2s2n'S2(n-l) “ 'S3n'S5(n-l) “ 'S5nS3(n-l) + S4n'S6(n-l) + S6n'S4(n-l)); 

Bn =-^(2sinS2(n-l) + 2S2n'Sl(n-l) + 'S4n'S5(n-l) + 'S5n'S4(n-l) + S3n'S6(n-l) + S6n'S3(n-l)); 

Bn 'S3n,S5(^^4) 'S5n,S3(^^4) + 'Sgn,S4(^^4) S4n'Sg(^^4), 

Bn —S4nS5(n+l) ~ '®5n'S4(n+l) “ 'S6n'S3(n+l) “ 'S3nS6(n+l)i 
Pn =2(sinS3(n+l) “ 'S3nSl(n+l) + 'S4nS2(n+l) “ 'S2nS4(n+l))! 

Hn =2(s2nS3(n+l) “ 'S3nS2(n+l) “ 'S4nSl(n+l) + 'SlnS4(n+l)) 5 
In '®5}^S3(^_4^ -j- S4nSQ(^n—l) '®6n'S4(n—1) > 

Jn 'S5^S4(^_4) 'S4^S5(^_4) 'S3,^Sg(^_4) + Sg^S3(_^_4), 

Bn 2(s3^S4(jj_ 4) '5ln'S3(n—1) “1“ ^2n-S4:{n—l) '54n'S2(n—1)) > 

Pn 2(s3^S2(n—1) '52n'S3(n—1) '5ln'S4(ri—1) “1“ '54n'Sl(n—1)) > 

Pn 2(s5^S4('jj_I_ 4) 'SlnS5(n+l) ^6nS2{n+l) “I" '52nSg(n+l)) > 

IPn 2(s5^S2(n+l) '52nS5(n+l) “1“ '56nSl(n+l) '5lnSg(n+l)) > 

Bn 2(s4^S5('jj_ 4) 'S5nSl(n—1) '52nSg(n—1) “1“ '56nS2(n—1)) > 

Bn =2(s2nS5(n—1) ~ '55nS2(n—1) “I" '5lnSg(n—1) ~ '56nSl(n—1)) • 

It is interesting to note that if set = fn,g-n — 9n /n = On + ibn,gn = Cn P idn, then we can 
show that the matrix Sn possesses the form 

•Sin i{s3n “ S2n) 

i{s3n + S2n) —Sin 

where the vector S„=(sin, .S 2 n, S 3 n)'^ G in i.e., + s^n “ •sin = Sjn,j = 1,2,3 is 

given by 

_ O^ndn bnCn _ Iffnl |/n| _ \fn\ Iffnl 

Q'nCn “1“ bndn ‘I(flnC-n “1“ bndn) 2(a,nCn “1“ bndn) 
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Thus, equation ([67)1 leads to 


<2 _o SnXSyt—1 ^ Sn,+lXSn, 

1-S„+1-S„- ^ ^ 

This means that the equation (I67p is a discrete modified Heisenberg equation under the condition 
/7„ = fn 5 9-n — 9n which agrees with a well-known fact that the discrete defocusing NTS is gauge 
equivalent to the discrete modified Heisenberg equation m- 

4 Conclusion and Discussion 

In this paper, we have shown that, under the gauge transformations, the nonlocal focusing NLS (the 
nonlocal defocusing NLS) and its discrete version, discrete nonlocal focusing NLS (discrete nonlocal de¬ 
focusing NLS) are, respectively, gauge equivalent to a Heisenberg-like equation (modified Heisenberg- 
like equation) and a discrete Heisenberg-like equation (discrete modified Heisenberg-like equation). 
From the gauge equivalence, we can see that the properties between the nonlocal NLS and its discrete 
version and NLS and discrete NLS have great differences. We have also obtained the discrete soliton 
solutions for the discrete nonlocal NLS through constructing the Darboux transformation. These dis¬ 
crete soliton solutions are different from ones obtained by using the scattering transformation. We 
should point out here that geometric interpretation for the nonlocal NLS and its discrete version is 
not very clear at the moment. How to understand the meaning of equations (1141) . ()25I) . (I48I ), and (|69p ? 
This problem is worth a further investigation in the future. 
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